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Abstract 

Contributions through second order, O(a^), in perturbative quantum chro- 
modynamics are calculated analytically for inclusive associated production 
of a prompt photon and a charm quark at large values of transverse mo- 
mentum in high energy hadron-hadron collisions. Seven partonic subpro- 
cesses contribute at order a^. We find important corrections to the lowest 
order, 0{as), subprocess eg 7c. We demonstrate to what extent data from 
p + p J + c + X may serve to measure the charm quark density in the 
nucleon. 
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I. INTRODUCTION 



Because photons couple in point-like fashion to quarks, observation, among the final- 
state particles in a high energy collision, of photons carrying large values of transverse 
momentum provides an incisive probe of the short distance hadron dynamics of the collision. 
This fact explains the substantial theoretical and experimental interest shown in studies of 
the cross section for production of photons at large angles in hadron-hadron and lepton- 
hadron scattering and in electron-positron annihilation processes. At stake are precise tests 
of the theory of perturbative quantum chromo dynamics (QCD) and use of data to determine 
properties of the relativistic proton such as the momentum distribution of its constituent 
gluons and quarks. Discovery of the charm quark and, later, of the bottom quark stimulated 
interest in the dynamics of their relatively copious production in high energy interactions of 
hadrons. Recent experimental advances now offer the possibility of studies of the associated 
production of a photon (7) carrying large transverse momentum along with a heavy quark 
(Q) whose transverse momentum balances a substantial portion of that of the photon. 
In this paper, we report a fully analytic next-to-leading order QCD calculation of the two- 
particle inclusive distribution for prompt photon plus associated heavy flavor production at 
large values of transverse momentum, with speciflcation of the momentum variables of both 
the flnal prompt photon and the flnal heavy quark. These results should facilitate further 
experimental tests of correlations inherent in the QCD matrix elements and provide a means 
for measuring the charm quark density in the nucleon. 

Although a qualitative description may be obtained from lowest-order perturbation the- 
ory, more precise predictions of the momentum distribution for the inclusive production a 
heavy quark (or antiquark) require perturbative calculations that extend to higher order. p| 
Likewise, perturbative QCD calculations of inclusive and isolated prompt single photon pro- 
duction are available. At the level of two-particle inclusive flnal states, next-to-leading 
order QCD calculations have been done for 77 production [|^,|^ , for 7-hadron production 
and for QQ correlations. ^ The cross section for the production of two hadronic jets has 
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been studied at 0{a^) by several authors. |T0| Constraints on the charm and strange quark 



densities from data on intermediate vector-boson production are discussed in Ref. ||11|| . 

For values of transverse momentum of the heavy quark significantly larger than the 
mass rriQ of the heavy quark, the cross section for the two-particle inclusive reaction p + 
p ^ 'y + Q + X may be calculated from the leading order QCD subprocess, the quark- 
gluon Compton process, g + Q 'y + Q- This subprocess is of first order in the strong 
coupling strength a^. The cross section is obtained as a convolution of the hard-scattering 
QCD matrix with probability distributions that specify the initial gluon and heavy quark 
constituent momentum densities in the incident hadrons, p and p. At next-to-leading order 
in QCD, several subprocesses contribute to the 7 + Q final state: 
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For computation of the cross section for Q production, the set of next-to-leading order 
subprocesses is obtained from those of Eq. (|1.1| ) after replacement of the initial Q's by Q's 
in Eqs. ( |l.la| ), ( |l.ld|) , ( p..le| ), ( p..lg| ). We note that for values of p^ that are comparable 



to or less than mq there would be no 0(0;^) subprocess, and the proper hard scattering 
expansion would entail only the subprocesses of Eqs. (|l.lb| ) and (|l.lc|) . For the remainder 
of this paper, we limit ourselves to charm production, and we work with the massless Q 
approximation, rric = 0. 

We are interested ultimately in the fully differential two-particle inclusive cross section, 
E^Egda/d^p^d^PQ, where {E,p) represents the four- vector momentum of the 7 or Q. For 
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each subprocess listed in Eq. ( |1 . 1| ) , this calculation requires integration of the momentum of 
the unobserved final parton {g, Q,q, or q) and over the initial parton momentum densities. 
CoUinear singularities are handled analytically by dimensional regularization and absorbed 
into initial-state parton momentum densities or final-state fragmentation functions. To make 
the analytic calculation tractable, we chose to work in terms of the transverse momentum 
of the final 7, p^, and the ratio of the heavy quark and photon transverse momenta: 



Q 7 
Pt-Pt 



;i.2) 



iPir ■ 

To warrant use of perturbation theory (and the massless Q approximation), we limit our 
considerations to z > and pj- > 10 GeV. The results should be applicable quantitatively for 

^ rric. The distribution in z from the leading order subprocess g + Q—^^y + Q is peaked 
sharply at 2; = 1 (a 6{1 — z) function in the naive collinear initial parton approximation). 
The next-to-leading order processes alter the size of this sharp peak and produce a broad 
distribution above and below z = 1. 

Contributions to hard photon production from long-distance quark to photon and gluon 
to photon fragmentation processes have been emphasized theoretically, |12| parametrized 
phenomenologically in leading order, |TE[ and evolved in next-to-leading order. |T5;|T3 These 



terms may account for more than half of the calculated inclusive single photon cross section 
at modest values of transverse momentum at the Fermilab Tevatron collider. Because of 
our kinematic restriction z > 0, there will be no contributions to the final cross section 
from (5^7 fragmentation, where Q is the observed quark/anti-quark, from among the 
subprocesses in Eq. ( |1 . 1|) . On the other hand, fragmentation of the unobserved final par- 
ton into a photon in subprocesses (l.a-g) will contribute to the cross section and produce 
photons that carry p^ less than that of p^, mostly populating the region z > 1. Photons 
originating through fragmentation are likely to emerge in the neighborhood of associated 
hadrons. An experimental isolation restriction is needed before a clean identification can 
be made of the photon and a measurement made of its momentum. Isolation reduces the 
size of the observed fragmentation contribution. To represent the effects of isolation, we 



should use fragmentation functions defined with a cone size. Photon isolation complicates 
the theoretical interpretation of results, however, since it threatens to upset the cancellation 
of infra-red divergences in perturbation theory. ^ In this paper, we calculate the contri- 
butions from photon fragmentation at leading order only, and, except for one illustrative 
figure, we neglect the isolation requirements. 

After integration over the longitudinal momentum of the heavy quark, we present our 
results in terms of the cross section da / dpj'dy^dz. Here, y'^ represents the rapidity of the 
7. Our desire to perform a fully analytic calculation restricts our ability to provide a more 
differential cross section in this paper (i.e., a cross section also differential in y^). In a later 
more detailed paper, we will present such results obtained from a versatile combination of 
analytic and Monte Carlo techniques. [0 In that method, selections may be made on several 



variables and photon isolation restrictions are easier to impose. An earlier theoretical paper 
addresses prompt photon plus associated charm production at large values of transverse 



momentum, as we do here, but our analysis differs from that of Ref. ||T^. The calculation 
of the photon-plus-charm cross section in Ref. is done in lowest order while ours is 



done at next-to-leading order. In lowest order, the subprocesses gg — > 7CC and qq — > 7CC 
contribute in the massive case, whereas eg 7c plus fragmentation processes contribute in 
the massless case. In a forthcoming paper, we intend to examine the massive case in detail 
and to discuss comparisons with the massless case in the regions of phase space of their 
respective apphcability. As remarked above, our massless approach should be appropriate 
and applicable in the domain in which there is effectively only one large scale, ^ m^. 

For the interval in pj. of current experimental interest, 10 GeV < < 50 GeV, the 
gc and gg subprocesses of Eqs. ( |l.la|) and ( |l.lb| ) are the most important quantitatively at 
Fermilab Tevatron energies, owing to the strength of the gluon density. For p7p > 70 GeV, 
calculations of the inclusive yield of single photons indicate that the qq subprocess begins 
to dominate, but the cross section is small in this region. Dominance of the perturbative 
subprocess initiated by gc scattering is preserved after the next-to-leading terms are included, 
justifying use of data from p + p^'y + c + X in attempts to measure the charm quark 



momentum density in the nucleon. However, we show that other subprocesses account for 
about 50% of the cross section at currently accessible values of pj^. The "background" 
associated with these subprocesses must be taken into account in analyses done to extract 
the charm density. 

Our results are provided in terms of the momentum of the charm quark. In a typical 
experiment, ^ the momentum of the quark may be inferred from the momentum of prompt 
lepton decay products or the momentum of charm mesons, such as D*'s. Alternatively, 
our distributions in z or may be convoluted with charm quark fragmentation functions, 
deduced from, e.g., e+e^ annihilation data, to provide distributions for the prompt leptons 
or L)*'s. 

In Sec. II, we present our analysis of the leading and next-to-leading order contributions 
to the partonic hard-scattering cross sections. Numerical results are described in Sec. Ill, 
and a summary of our conclusions is provided in Sec. IV. An Appendix is included in which 
we present our method for performing the required three-particle final-state integrals in 
n-dimensions to extract the singularities of the two-particle inclusive hard cross section. 



We consider the two particle inclusive reaction A + B ^y + c + X where A and B denote 
incident hadrons; p'^ and denote the four- vector momenta of the photon and charm quark. 
The usual Mandelstam invariants are defined in terms of the momenta of the two incoming 
hadrons P4 and Pb, and the momentum fractions of the initial partons, Xi and X2, via 



II. ANALYTICAL CALCULATION 
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Here ^/s is the center-of-mass energy in the hadronic system. We define 



V = l + - 
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A. Leading Order Contributions 

In leading order in perturbative QCD, only one direct subprocess contributes to the 
hard-scattering cross section, the QCD Compton process eg —>■ 7c, unlike the case for 
single inclusive prompt photon production, where the annihilation process qq — > 'jg also 
contributes. Since the leading order direct partonic subprocess has a two-body final state, 
the photon and c-quark are produced with balancing transverse momenta, and the variable 
z, defined in Eq. (|1.2|), is always unity. 

The leading order direct partonic cross section is 

"^"^ ^6{l-z)5{l-w), (2.3) 



dvdzdw dv 
where da/dv is the partonic Born cross section: 

da 1 vraema^e^ 1 + (1 - t;)2 

— {eg ^ 7C = ^ 2.4 

dv Nc s 1 — 1> 

Here and are the electromagnetic and strong coupling constants, respectively, Nc = 3 
is the number of colors, and Cg denotes the quark charge. 

The full expression for the physical cross section in leading order is 

- ^"""^i L T^,f'^^- ""^f^^^^- "'^f/^'- ^'^"(^ - ■"■) + P-^' 

Quantities V and W are defined similarly to v and w, Eq. (|2.2| ), but in the hadronic system; 
f^{xi, M^) denotes the parton density in hadron A as a function of the momentum fraction 
Xi and factorization scale M. 

In addition to the lowest order direct subprocess just discussed, eg 7c, there are 
fragmentation contributions that are also effectively of leading order in a.^- In these contri- 
butions the photon is produced through fragmentation of a final-state parton from any of 
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the subprocesses listed below. The fragmentation functions are essentially of order 

0{aem/as) 
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(2.6) 

We are interested in configurations in which the photon and charm quark have relatively 
large and to-some-extent balancing values of transverse momentum. Therefore, in the cases 
of the first, third, and fourth of the subprocesses listed above, the photon is produced from 
fragmentation of the g and non-charm quark q, respectively. In the other cases it is produced 
in the fragmentation of one of the (anti) charm quarks. The expression we use to evaluate 
the fragmentation contributions is 

da . 1 dz' /•! dv ,.o..„n, .,2n^^^"^"*^ 



apj^ay^az its Ji-v+vw z'^ Jvw I — v av 



xD,/,iz',Q')5i^-z'). (2.7) 



In a fully consistent next-to-leading calculation, one should calculate the subprocesses in 
Eq. ( |2.6|) to 0(af), since the photon fragmentation functions that are convoluted with the 
hard subprocess cross sections are of O^aem/ds)- For simplicity, we include them in 0{al) 
only. In fact, next-to- leading order fragmentation contributions to single prompt photon 
production have been included only once before [Q. We expect the next-to- leading order 
corrections to the fragmentation contributions to be insignificant numerically especially after 
isolation cuts are imposed. 
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B. Next-to-leading order Contributions 



There are two classes of contributions in next-to-leading order. First there are the virtual 
gluon exchange corrections to the lowest order process. Examples are shown in Fig. lb. These 
amplitudes interfere with the Born amplitudes and contribute at O(aemtts)- They have been 
calculated twice before. We use the results of Ref. The virtual contributions are 
proportional to 6{1 — z) and 6{1 — w). At next-to-leading order there are also three-body 
final-state contributions, listed in Eq. (|1 . 1| ) . The matrix elements for these are also taken 
from Ref. ||^, where they are calculated for single inclusive prompt photon production. 

The main task of our calculation is to integrate the three-body matrix elements over the 
phase space of the unobserved particle in the final state. The situation here is different from 
the standard case of single inclusive particle production, first developed in Ref. [l^, since we 
wish to retain as much control as possible over the kinematic variables of a second particle in 
the final state, while at the same time integrating over enough of the phase space to ensure 
cancellation of all infrared and collinear divergences, inherent when massless particles are 
assumed. Because our goal is to provide a fully analytic calculation, we find it necessary 
to integrate over the full range of rapidity of one of the observed final-state particles. We 
choose to integrate over that of the charm quark, since the photon is usually considered the 
trigger particle in the experiments. 

The situation here is similar to that met by Aurenche et al @,|], and we use a similar 
technique to perform the phase space integrals. We give a fairly detailed outline of the 
method since it is necessary to adapt it to our situation and also because it has not been 
widely used. We believe our presentation clarifies certain details which are not stressed in 
the above references. 

The three-body phase space integration is done in the rest frame of the observed c 
(or c-quark) and the third unobserved parton. Denoting the momenta of the process by 
Pi + P2 — ^ ^1 + ^2 + ^3, we work in the rest frame of k2 and fcs, where ki is the momentum 
of the trigger photon. The final form of the three-particle phase space integral (see the 
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Appendix) is 
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(2.8) 



1 — vw 



We are left to perform the final integration of the squared matrix elements over 62. 

As in the case of single inclusive cross section calculations, documented extensively else- 
where, one can use relations among the Mandelstam variables to reduce complex combina- 
tions of them to simple products and ratios. The phase space integral over 62 is performed 
in 4 — 2e dimensions, thereby exposing coUinear and soft singularities as poles in e. After 
the three-particle phase space integrals are performed, we obtain a three-body final state 
hard-scattering cross section that we represent by the expression 



Superscript R indicates that this is the subprocess cross section for a real three-body final- 
state contribution, as distinct from the contribution from the virtual gluon exchange con- 



In general, has single and double poles in e. In accord with the factorization theorem 
of perturbative QCD, the double and some of the single poles cancel between the real and 
virtual contributions. The remaining single poles in e represent collinear divergences that 
are subtracted into parton densities and fragmentation functions. 

In order to illustrate how the collinear singularities are handled we discuss a few repre- 
sentative examples. 

(a) c + g —> + c + X 

This is the QCD Compton process plus higher order corrections. We label the momenta 




tributions that we denote aYj- The subscripts ij designate one of the processes in Eq. (p..ip. 



by 




(2.9) 
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In performing the phase space integration, we expect to encounter singularities where the 
gluon ^3 becomes soft and/or parallel to pi,P2 or k2- Since we require that the observed 
charm quark and 7 be in opposite hemispheres, we will not encounter any singularity where 
ki and k2 are coUinear (see the Appendix). In the cases where the gluon is either soft and/or 
parallel to pi or p2, then z — 1. We expose the z ^ 1 singularities by using the expansion 



-S{l-z) + -^ -^+ / .^ -2e e{l-z) + 0(e'). 



|l_^|i+26 _2e ^ ^ (1-^)+ (^-1)+ V 1-^) 

(2.10) 

There are plus-distributions in the variable z, as well as the usual ones in w that arise in the 
single particle inclusive case and correspond to the gluon becoming either soft or coUinear 
to ^2- Plus-distributions in z and w can be encountered simultaneously and must be treated 
carefully in the numerical evaluation of the cross section. 

Once the phase space integrals are performed and the soft and coUinear poles are exposed, 
we can add the real three-body contributions to the virtual gluon exchange terms, after which 
all the double poles cancel along with some single poles. The remaining coUinear poles must 
be factored into the parton distribution and fragmentation functions. We perform these 



subtractions in the universal or MS scheme, described in detail in many places. 

To account for all coUinear configurations allowed in the subprocess, the counter cross 
section or factorization formula that must be added to our results in order to cancel the 
coUinear poles is 



1 d<T^ a. 
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sv dv 



+ ^^99 , — 1 e)5(l - Z) 

s[l — vw) \1 — vw ) dv 

+ ^H,c{z, M"')^ (s, V, e)e{l - z)5{l - w) 

sv dv 



(2.11) 



H,j{z,Q^)^-jPij{^) 



+ fij{z), (2.12) 



and 
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(2-13) 

Functions Pij{z) are the one-loop splitting functions [18], fij{z) = and dij{z) = in 



the MS factorization scheme, and fx is the renormalization scale. In the MS scheme, 
1/e = l/e-'jE + InAn. 

In Eq. ( p.ll|) , we distinguish the factorization scale M and the quark to quark plus gluon 



fragmentation scale M". The last term indicates that we factor the coUinear singularity 
that arises when the observed charm quark k2 becomes parallel to the gluon, k^, into a 
fragmentation function at scale M"^, for the production of a charm quark. Note that this 
singularity occurs in the region z < 1, since the photon must balance the momentum of the 
charm-gluon system. 

We are free to convolute our cross section with a fragmentation function that describes 
the formation of specific charm decay products (e.g., D or D* mesons), but we choose not 
to do so in this paper. 

(^)9 + 9-^l + c + c 

In the gluon-gluon fusion process, gg — 7cc, the photon may become coUinear to the 
unobserved final-state quark, a situation not encountered in the gc process discussed above. 
This singularity occurs a.t z = zi where zi = — v + vw), and, as discussed in the 
Appendix, we use an expansion similar to that in Eq. ( p. 10 ) to expose the singularity. Note 



that this singularity occurs in the region z > 1, and that z is exactly the reciprocal of the 
usual fragmentation variable for a parton to fragment into a particle with a fraction of its 
momentum, 1/z. The factorization formula for this process is 



1 da^ 



a. 



sv dvdwdz 27r 



1 da'^^^^'^ 

—H,g{w, M^)— {ws, V, e)6{l - z) 

sv dv 



+ TT^ ,H.J-, A iws,vw,e)5il-z) (2.14) 

s{l-vw) \1 - vw J dv 

+ \ 1 -v + vw, M'2 — s, ■ , e)S{zi - z) 

s[l — v + vw) dv 1 — v + vw 

In this equation, we distinguish the factorization scale M and the quark to photon fragmen- 
tation scale M' . 
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(c) q + q^j + c + c 

The process qq —>■ '-fcc, as well as that of Eq. ( |l.ltD , has a final-state collinear singularity 
when a gluon splits into a collinear cc pair, and, in addition, a singularity when the photon 
is produced from fragmentation of a final-state quark. The factorization formula for this 



case IS 



sv dvdwdz 



2^ 



sv 



dv 



+ 



1 



s(l — f + vw) '^^ 
X 6{zi - z)] . 



H^sil -v + vw, M'^) 



da 



qq—tcc 



VW 



dv 



V + VW 



(2.15) 



C. Physical cross section 

Once all singularities are dealt with, we calculate the physical cross section by convoluting 
the hard partonic cross section with parton distribution functions. In terms of the variables 
we are using, the cross section at next-to-leading order is 

da 



.1 ^o.-y^^ ^ jyy^ Y — V JvW/v W ■' 

5(1 - z)5{l -w) + ^^^^^Kij{s, V, w, z, fi^, M^, M'\ M"^) 



dpji'dy^dz TTs ~^ Jvw 1 — v Jvw/ 



V dv 271 
The first term within the square brackets is the leading order part, and 



(2.16) 



Kij{s, V, w, z, /i', M\ M'\ M'") 

is the next-to-leading order correction term; K^j may include virtual gluon exchange contri- 
butions. 

Taking the eg subprocess as an example, we outline how we obtain the function 
Kij{s,v,w,z,fj,'^,M'^,M'^,M"'^). The virtual gluon exchange contributions are represented 

by 

2 1 1 



dvdwdz V e 
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They are proportional to S{1 — w) and S{1 — z). The real three-body contributions are 
denoted 

^< 1 1 

avawaz \ e 

Combining the three-body final-state contribution and the virtual gluon exchange contri- 
bution and adding to these the subtraction term in Eq. (|2.11| ), we derive a finite subprocess 
cross section: 

K,gis, V, w, z, /i^ M\ M"2) = ^^^^ (^s, V, /i^, 1 1) + 

' s,v,w,z,—,-] + 



dvdwdz \ e. 

-^Uv,w,z^,M\M"A. (2.17) 
dvdwdz \ e J 

At this stage all single and double poles cancel, and we are left with a finite cross section 
dependent on the factorization scale M and fragmentation scale M" . Because of the ad- 
ditional variable z, the function K^g is quite lengthy when compared to that for inclusive 
single photon production. In schematic notation, where only the 2;-distributions are made 
explicit, we can write the hard-scattering cross section as 

K,g{s, V, w, z, fi\ M\ M"2) = ci{v, w)5{l - z) + C2{v, w) f|^~ f 

[I - zj+ 

. Mz-V) . ,/ln(l-^)\ 
+ c-i[y,w)- — — V C4{v,wj 



[z-l)+ ' ' \ 1-z 

+ C5{v,w,z). (2.18) 

The functions Ci{v, w) contain, in general, distributions in (l — w), and they can be expressed 

by 

c,{v, w) = c]{v)S{l -w) + (^^^ + cKv) (^^^) 

+ ct{v,w). (2.19) 

Similar expressions can be written for the other subprocesses. These will generally involve 
the fragmentation scale on the photon leg, M', and additional distributions in {zi — z) and 
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{z — Zi). These are defined as normal plus-distributions, but in the intervals [0,zi], and 
[zi.z^ax] 5 respectively. We integrate the distributions between limits other than these. For 
example, if the limits in the first case are [za-, -Zi], we must make the replacement 

1 1 



{zi - z)+ {zi - z),^ 
where the new distribution is defined by 



+ 6{zi- z)\n{zi-Za), (2.20) 



''dz^I^= TdzM^I^. (2.21) 

[Zi — Z)z^ Jza Z\ — Z 

By expanding our integrated matrix elements as plus-distributions in z^ we are able to 
expose the singularities that occur at z = 1 and z = z\. This procedure ensures that these 
integrable singularities can be treated numerically. However, it also means that our analytic 
distributions in z are singular at z = 1 and z = Zx. For comparison with experiment, we 
provide predictions for the z dependence in the form of histograms with finite bin-widths 
Az, reminiscent of experimental resolution. As in Ref. we define 



da 1 da 



dpjpdy^dz Az Jz-^ dpj^dy^dz' 
For distributions in pj,, we integrate over a specified range of z 

da r^b da 



dz'. (2.22) 



dz. (2.23) 



dpj^dy^ Jza dpj^dy^dz 

This completes our discussion of the calculation. Further details can be found in the 
Appendix. 



III. NUMERICAL RESULTS AND DISCUSSION 

In this section we present and discuss explicit evaluations of the correlated production 
cross section of charm plus a prompt photon. We provide results at pp center-of-mass energy 
y/s = 1.8 TeV appropriate for the CDF and DO experimental investigations underway at 
Fermilab. The cross sections we evaluate are those derived in the text: Eqs. (|2.5|), ( |2.7] ) and 
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( p.l6| ). For the electromagnetic coupling strength we use aem = 1/137, and we employ a two- 
loop expression for ^^(/i^) with quark threshold effects handled properly. We choose identical 
values for the renormalization, factorization, and fragmentation scales, fi = M = M' = M" . 
In the results presented below, we vary /i to examine the sensitivity of the cross section to 
its choice. We choose Ag^^ according to the parton distribution set we use; J^^qco = 0.200 
for the GRV parton distributions. The sums run over 4 flavors of quarks (m, (i, c, s), all 
assumed massless. We do not include a h quark contribution in our calculation. 

Most of the calculations reported here are done with the GRV parton densities . We 



observe some differences when we use instead the CTEQ3M densities [^. The magnitude 
and Bjorken x dependence of the charm quark density in these two sets are similar, as 
shown in Fig. 2, but show some differences at large x, leading to a 30% difference in the 
cross section at pj^ = 60GeV. In these densities, the charm quark probability is generated 



com- 



through perturbative evolution, and there is no non-perturbative intrinsic charm |2^ 
ponent. Neither density may be correct since there is little direct experimental information 
to constrain this density [|ri[|. A goal of our analysis is to ascertain the extent to which the 
gc initial state is expected to dominate the cross section for p + p — > 7 + c + X, and, thus, 
the extent to which data from this reaction may serve to measure the charm quark density. 
The quark-to-photon fragmentation function is expressed as 



2x "e 



.2.21-1.28^ + 1.29^2 ,2 1 54 

Co ^ ^ 2 + 0.002 1 - z? z~^-^^ 

^ l-1.63£n 1-2 ^ ' 
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xin^fi'/fil). (3.1) 



The gluon-to-photon fragmentation function is 



z Dg^Jz, /i^) = ^ 0.0243 (1 - z) z-^-^^ in I [it] • (3.2) 
27r ^ ' 



These expressions for -D„^^ and D„^^, taken from Ref. |13|, are used as a guideline for 



our estimates. The physical significance of scale /io is that the fragmentation function 



vanishes for energies less than /iq. For the m, c?, s, and c quarks, we set /io = Ag^^ 



as 



in Ref. ||T^. We remark that we use simple leading order fragmentation functions in our 
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calculation, in contrast to the fact that we have done a next-to-leading order MS calculation. 
It would be more consistent and, therefore, preferable to use MS fragmentation functions 
evolved in next-to-leading order. Our choice of leading-order fragmentation functions is 
motivated by our desire to work with analytic expressions. In published analyses of next- 
to-leading order fragmentation functions, |]15|15| the general formalism is presented but the 
fragmentation functions themselves must be obtained through numerical evolution codes. 
Our primary purpose in this paper is to provide a theoretical framework for the analysis 
of the correlated production of charm and prompt photon not necessarily to present the 
most up-to-date numerical predictions. Thus, we believe our leading-order fragmentation 
functions are adequate. 

In several figures to follow, we show the predicted behavior of the photon yield as a 
function of and z, as well as the breakdown of the total yield into contributions from the 
leading order and the various next-to-leading order pieces. The ratio z is defined in Eq. ( [1.21) . 
We choose to display cross sections as a function of the ratio z, for fixed values of or as 
a function of pj^. We choose the renormalization/fragmentation scale = pip. Since both 
the photon and final charm particle carry large transverse momentum, we could perhaps 
equally well choose fi = p^ or some combination of the two. In selecting pj,, we focus upon 
the photon as the "trigger" particle whose transverse momentum is well determined. We 
display the /i dependence of our results below. 

Throughout this paper, for clarity and simplicity of the discussion, we refer consistently 
to charm production, e.g., p + p— >7 + c + X. However, the numerical values of the cross 
sections shown in the figures are those for the sum of charm and anticharm production in pp 
scattering. In Fig. 3, we present the photon yield as a function of the ratio z for two choices 
of pj'. The same results are displayed in Fig. 4 as a function of p^ for z integrated over the 
interval 0.2 to 2.0. We restrict z > 0.2 as otherwise the transverse momentum of the charm 
quark could become unacceptably small. In Fig. 3(a), the net lowest order contribution is 
shown at p^ = 15 GeV. The lowest order contribution is made up of the lowest order direct 
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term, eg 7c, and the fragmentation terms discussed in Sec. II. A. The direct term provides 
a 5-function at 2; = 1 since the photon and charm quark carry equal but opposite transverse 
momenta at this order. The parton to photon fragmentation contributions populate the 
region z > 1. In the collinear fragmentation, the photon's transverse momentum is opposite 
to that of the charm quark but its magnitude is less. One of the striking features of Fig. 3a, 
is that the net fragmentation contribution to the cross section is quite small compared to 
the case of inclusive photon production. At Tevatron energies, fragmentation accounts for 
about 50% of the inclusive yield at this value of pt [HHil- (Note that we have not yet 



imposed any isolation cuts on the cross section.) One reason for the small fragmentation 
contribution is that fragmentation from the eg initiated process is strongly suppressed due to 
our restriction that the charm quark and photon be in opposite hemispheres (-2 > 0). Thus 
only fragmentation from the gluon leg is included, and the — >■ 7 fragmentation function is 
in general smaller than that for 5^7. 

In Figs. 3(b) and Fig. 3(c), we show the z distribution after the next-to- leading order 
contributions are included. The solid lines show the full result in which both the lowest 
order and all next-to-leading order terms are incorporated. Comparing the solid curve in 
Fig. 3(b) with that in Fig. 3(a), we note that the z distribution is substantially altered once 
the next-to-leading order terms are included. In particular, the peak at z = 1 is reduced in 
magnitude by about a factor of 2, and the z distribution gains significant breadth below and 
above z = 1. The reduction in the magnitude of the peak at z = 1 is attributed to the effect 
of the O(a^) collinear contributions on the initial parton legs. These collinear terms provide 
the same event structure as the lowest order direct subprocess, viz., a final-state photon and 
charm quark with equal but opposite transverse momenta, but their contribution is negative 
due to ln(l — z) terms from the phase space and large logarithms of (1 — Zmin) and {z^ax — 1) 
from the 1/(1 — z)+ and (2; — 1)+ distributions; Zmin and z^ax are the lower and upper edges of 
the bins around z = 1. On the other hand, away from collinear configurations, the 0{aemCel) 
subprocesses, listed in Eq. (|1.1D, generate three body final states in which three final partons 
share the transverse momentum balance. The non-collinear contributions therefore populate 



a broad interval in z. 

In addition to the complete result through next-to-leading order, the solid line in Figs. 
3(b) and 3(c), we display also contributions from three of the 0{q?^ terms. The sum of the 
contributions from the other four 0{p?^ terms is negligible by comparison at = 15 GeV. 
The individual contributions show the important role that the O(a^) terms play at values 
of z below and above 1. Contrasting Figs. 3(b) and 3(c), we see that the peak near z = \ 
is predicted to sharpen as is increased, reflecting a diminishing importance of the 0{p?^ 
terms at larger transverse momentum. 

In Fig. 4, we show the cross section as a function of the transverse momentum of the 
photon, To obtain these results, we integrate over the interval 0.2 < z < 2.0. These 
results show that the eg intial state dominates the cross section until pj- approaches 100 GeV. 
It accounts for 60%, 55%, and 50% of total at pj' = 15, 45, and 60 GeV, respectively. The 
gg contribution is important at small values of pj-, but it falls off more steeply with p^i than 
the eg contribution. The contribution from the valence subprocess, qq ccr)^ is negligible 
at small pj,, but it overtakes the contribution of the eg subprocess at sufficiently large pj,. 
Owing to the fact the valence quarks carry significantly harder fractional momentum than 
the gluons and charm quarks, a major role for the valence subprocess is expected at large 
enough pj.. However, the numerical results indicate that the hard-scattering matrix element 
overcomes this effect at modest values of pj., resulting in dominance of the eg initial state. 
Comparison of Fig. (4) and Figs. 3(b) and 3(c) shows significant z variation in the fraction 
of the total cross section accounted for by various subprocesses. 

Dependence on the renormalization/factorization scale // is displayed in Figs. 5 and 6. 
As \x is increased, ag decreases, resulting in a reduction of the hard-scattering cross sections. 
The parton densities also steepen as \i is increased. Both effects contribute to the typical 
decrease of the cross section at fixed large pj- as is increased, as shown in Fig. 5. The 
\x dependence of the z distribution presented in Fig. 6 is considerably more significant. The 
distribution becomes more sharply peaked at 2; = 1 as is increased. As shown in Fig. 3(a), 
the leading order direct contribution produces a sharp peak at 2; = 1, whereas the next-to- 

19 



leading order contributions broaden the distribution, as shown in Figs. 3(b) and 3(c). The 
decrease of as n increases diminishes the relative importance of the next-to-leading order 
contributions. 

The functional form of Dg^^ {z, /i^), Eq. ( |3.1| ), shows that the fragmentation contribution 
increases logarithmically as /i is increased. If the fragmentation contributions played a major 
role in the final answer, one would expect different /i dependence from that shown in Fig. 6. 

In Fig. 7, we present the "i^-factor" as a function of py. Here K is defined as the 
ratio of the complete answer through next-to-leading order to the full leading order answer 
(including the leading order fragmentation terms). Our results show that for z > 0.2, 
the inclusive K factor is about 2 for pj. > 15 GeV. In the inclusive case, no isolation 
requirement is imposed on the photon. To make contact with experiment, an isolation 
restriction is necessary. Because fragmentation contributions do not play a significant role 
in the associated production of photon plus charm for z > 0.2, we do not expect a great 
change of the K-factor after isolation is imposed. To estimate the impact of isolation, we 
use a combination of analytic and Monte Carlo methods. |T^ We choose an isolation cone 
size R = 0.7, and energy resolution parameter, e = 2 GeV/pj., as is done in the CDF 
experiment [|I|. We find that the K-factor is reduced to about 1.5, in respectable agreement 
with experimental indications. ^ 

IV. CONCLUSIONS 

In summary, we have computed the contributions through O(a^) in perturbative QCD 
for inclusive associated production of a prompt photon and a charm quark at large values of 
transverse momentum in high energy hadron-hadron collisions. The next-to-leading order 
terms alter the expected distribution in the ratio of the magnitude of the transverse momenta 
of the charm quark and prompt photon in an interesting and measurable fashion. The overall 
cross section increases by about a factor of two after the next-to-leading terms are included. 
Dominance of the perturbative subprocess initiated by gc scattering is preserved after the 
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next-to-leading terms are included, justifying use of data from j9-|-p — > 7-l-c-l-X in attempts 
to measure the charm quark momentum density in the nucleonHowever, other subprocesses 
are shown to account for about 50% of the cross section at currently accessible values of pj-i 
and the "background" associated with some of these subprocesses, which are not initiated 
by charm quark scattering, such as in Eqs. ( |1 . 1 b| ) and ( |l.lc| ) must be taken into account in 
analyses done to extract the charm density. 
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APPENDIX A: THREE BODY CROSS SECTIONS 



In this Appendix we present a fairly detailed description of the techniques for performing 
the 3-body phase space integrals in n dimensions. We label the momenta for the general 
process by pi +P2 ki + k2 + k^, where pi and p2 are the incoming partons, and ki and 
k2 always label the observed photon and charm quark respectively. We integrate over the 
kinematic variables of k^. 

The calculation is performed in the rest frame of k2 and k^. In this frame of reference 
h + ks^O . 

Sij — {ki -\- kj) 
U = (Pi -Pif 

Ui^{p2-Pif, (Al) 

where i,j — 1,2, 3, and ti — t and ui — u, as defined in Sec. II. In terms of the momenta, 
the variable z is 

~- "'-"^ -m.k2, (A2) 



where m is a vector that depends on the choice of axes. We choose our axes in n dimensions 
such that 

/ s 

m = W— (sinhx, 0, 0, coshx) . (A3) 
V tu 

The axes are fixed and cannot be changed to simplify any phase space integrals we may 
encounter because, unlike the case of single inclusive particle production, we will not inte- 
grate over the full range of angles. The momenta of the particles can be parametrized in 
this frame as 

sv 

Pi = 7r-^(l,0, ...,0,sinV^',cos^') 

^^yS23 



s(l — vw) 
P2 = — (1,0, ...,O,smV;,cos-0) 

^^yS23 
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!(l^p^(l,0,...,0,rinr,co=r) 

^^^^^1, 0, 0, 806*1 008 6*2, 0086*1) 

1, 0, 0, — 8in6'i 008 6*2, — 0086*1). 



2_ 

S23 



ki 
k2 

Quantitie8 v and w are defined in Seo. II. 

From the definition of m, Eq. ( |A3D , we oan derive the relationships 



lw{l — v) 



tanh X = w 1 

V 1 — vw 

cos ip = cos ip' = tanh x 



sin ip = — sin ip' 



1 — w 
1 — vw 



cos ip" 
sin ip" 



1 + f — vw 
1 — V + vw 



tanhx 



1 — f — vw 1 — w 



1 — V + vw \ 1 — vw 



The oonstrained three-partiole phase spaoe is expressed as 

x5+{kl)5+{kl)5+{kl)5 (^-l-^-^ 

x5 \w -\ I 5{z — m.k2). 

\ s + t 



(A4) 



(A5) 



(A6) 



After some of the integrals are done with the aid of the (5-funotions, the element of phase 
spaoe reduoes in = 4 — 2e dimensions to 



P5(3) 



7TS /An 



Ait 



(27r)5 \ s J r(l - 2e) \swv{l - v) 
1 — w + Aw(l — v)z(l — z 



v-'(l-w)~'2, 



w{l — v) 



\(1 



vw 



X 



1 — vw 



(102 sin-2^( 



(AT) 



In partioular, we integrated over angle 61 using the funotion 6{z — m.k2) and the relation 



z = -{1 — 00s 6*1 00th x). 



(A8) 
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that can be derived from it. We are left with the task of integrating the squared matrix- 
elements over angle 92- 

Using relations among the Mandelstam invariants, and partial fractioning, we reduce 
functions involving 62 to only a few types for all subprocesses of interest. We denote the 
general invariant by 

Ti = Tio{ai + /3i cos 6*2), (A9) 

where ctj and f3i are functions of ■?/', and ip" , and hence of v and w\ TiQ is also a function 

of the latter (see Eqs. (|A4|) and ([A5|)). The combinations we must consider are 

IT" ^1 
— , 7^, and 



where, i,j,n = 1,2,3. These, in turn, are all expressible in terms of two general integrals, 
but the form of the functions a and (3 determines the final result, such as its singularity 
structure. 

The two general integrals are 

/o[T,] = h= r sin-2^ 92(192 = 7r2^^ ^i!r ^1 , (AlO) 
JO i ^[1 — ej 



and 



/i[T,] =T,o / 
Jo 



^sm-'^'92d92 r sm-^'92d92 



TT 



Ti Jo {a + (3 cos 6*2) 



a2 -/52 



r2[l -e] V2 ' ' a 



In terms of Jq, the following powers and combinations of propagators yield 

T^ =^ lo, 
Ti =^ TiQ^ailo), 

TiTj =^ TioTjo (^aiaj + ^^y^^) ^O' 
T! ^ Tl [a^ + h. (A12) 
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In terms of Ii we obtain: 



1 ^ 1 



T 

rp2 
3 



TT 
j-iJ-j 



iO 



Wo \ 



_(3, 



Pi 



rp2 



iO 



3 _ 



1 





2 

h[Ti]+7r 


[ A J 




3 

h[Ti\ +7r 


[ A J 



Wf 



1 



-imn-PMTj]). 



(A13) 



TioTjo aj(3i - ai(3j 

In order to demonstrate how the different propagators are handled in the calculation, we 
consider a few typical examples. We examine single propagators first, then double propaga- 
tors. 



1. Single Propagators 

In this case a — 1 — cos cos ^i, and (5 — —sva.%1) sin 9i. 

q;2 - ^ (cos Oi - cos ^^)f = Az'^ tanh^ x- (A14) 

There is a singularity at 2; = 0, but the physical condition that the photon and charm quark 
be in opposite hemispheres guarantees z > Q. The integral is therefore finite and can be 
treated in 4-dimensions. The result is 

Jo U2 2u2o tanh x z 

The result for t2 is similar. In the evaluation of U2 in terms of angles, 1*20 is the overall factor 
that does not depend on angles. 

(ii) ^ 
S12 
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Here a = 1 — cos ■i/'" cos ^i, and (3 = — sin -i/;" sin . 



(A16) 



1 — V + vw ^ 

A singularity occurs only for negative z when the photon and charm quark are exactly 
collinear and in the same hemisphere. We can treat this integral in 4-dimensions for the 
same physical reason as above, with the result 

r^sm-^'92d92 vr 1 1 



/ 

Jo 



Si2 



2si2otanhxf2 + -H(i-^V 



(A17) 



[ill 



Ms 



In this case a = 1 + cos cos 6*1, and /? = sin sin 6*1. 



(A18) 



There is a singularity when z = 1, corresponding to and p2 being collinear. This pole 
must be exposed and factored into the parton distributions. The integral is 



/ 



^ sin"^^ 92(192 



IT 



U3 



1 + tanh^ x(l - 2z) 



2M3otanhx |1 — ^r^+^e 



tanhx 



r(i-26) 

P(l-e) 



(A19) 



This integral can be reduced if we use the expansion 
1 _ 1 ^(1 _ ,) + »(1 - . 



11 



-2e 



1-^ 



(A20) 



and note that the hypergeometric function at 2; = 1 reduces to 2 We introduce a plus- 
distribution, 



^max 

1 h-i] 



-dz 



Z^a. f(^z) - /(I) 



dz. 



(A21) 



where Zmax = 1/2(1 + cothx) from Eq. ( |A8D . When the phase space factor involving z is 
included in the expansion, the integral reduces to 
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sm 



62(162 



r(l-2e) 



Ms 



2^30 tanh x r2(l — e) 



5{l-z) 



In z. 



, 6{l-z) , g(^-l) 

max 11/^ \ 1 



{1-zU (z-l)- 



(A22) 



The result for l/ts is similar. 



Sl3 



This propagator occurs when there is a quark in the final state that may become collinear 
with the photon ki, such as in Eqs. ( |l.lb|) - (|l.lg| ). Here, a = 1 + cos cos 6*1 , and /? = 
sinip" sin 61. 

1 



a" - = 4 hz - 



tanh X- 



(A23) 



1 — f + , 

There is a singularity at z = 1/(1 — v + vw) = Zi. Using an expansion similar to that in 
Eq. ( |A20|) , but with {I — z) and {z — 1) replaced by (21 — z) and {z — zi), we cast the result 
in the form 

sin-^' 62(162 vr r(l - 2e) 



/•Tl 

Jo 



Sl3 



2si3otanhx r2(l - e) 
6{zi - z) 6{z - Zl) 



6iz-zi)[--^-\n (^1-^ 



(A24) 



{zi - z)+ {z - Zi)+_ 

Here Zmin = 1/2(1 — cothx), and the plus-distributions are defined in a similar way, but 
with limits from to Zi and Zi to Zmax- 

In this outline, we have omitted phase space factors present in Eq. ([A7|) and included 
them only when they are important for the expansions performed. In principle, all phase 
space factors should be included and, for example, v"^ would be expanded as 

e2 

1 - eln(f) + —\n'^{v) 

and combined with the final results of (i) to (iv) above, before e is set to zero. Factors such 



as 2v tanh x have also been omitted but are included in our final results. 
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2. Double Propagators 



We examine a few important examples of double propagators. Some have been calculated 
previously, |^,^ but there are cases not encountered in earlier calculations that we stress 
here. We include all phase space factors since most are needed in the expansions. 

1 



t3U3 



This double pole propagator was encountered in Ref. 0. The result is 



Jo 



sin"^^ 92d92 



nv 



r(l-2e) 



tgoUso r2(l - e) 



X 



(1 — vw) 



9{l-z) 



9{z-l) 



XD^il-z)U Pi(^-l))- 
/ 1-vw /tanhx + cothx\ 

+ ' [-T^ [ 2 j 

+ V (^-- - ln(l -vw) + 2 ln(l - w)^^ 



{l-v){{ - +ln(l 



'ln(l — w) 
(1-w) 



+ 



1 — vw 



1 — w 



In 



w 



6{1 - z)6{l - w) [ J-^ + 1 ln(l -v) + ^ ln\l - v] 



(A25) 



where Di = 1 — w + 2(1 — z)w{l 



[II] 



hU2 



Applying the last result in Eq. 



we get a term 1/(1 — 2z) multiplying the integrals for 



the propagators {[3{u2)I[u2\ + (3{t^)I\t?^), along with other factors. This term is singular at 
z = 1/2, but the singularity has no physical origin and must be removed before numerical 
evaluation of the cross section. The integral of l/u2 yields a term 1/z (Eq. (|A15|) ) while 
that of ts yields plus-distributions in {1 — z) and {z — 1) (Eq. ( |A22| )). To remove the false 
singularity we make the replacement 

1 ^(1 - z) 9{z-l) 

z z z 
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These 9 functions can be combined with the plus-distributions from the second term in 
{I3{u2)l[u2]+ I3{h)l[h]) to produce 



sm-'' 92(192 Tiv^-'{l-w)-'V(l-2e) 



hU2 



2^30^20 tanh^ X r2(l - e) 



5{l-z) 



In 2;„ 



+ 



'l-z) 



zl-z) 



+ 9(z-l] 



l + 2z 



>-!)+ (l-2-2)-2, 
The term 1/(1 — 2z) is harmless when multiplied by 9{z — 1). 

{iii) 



(A26) 



^3-523 



This case involves singularities when z ^ 1 and w ^ 1. It is discussed in Ref. The 
result is 



^sm-^'92d92 nv^-'T(l-2e 



^3^23 



+ 



svho r2(l - e) 

ln{z) f\n{l - z) 
^l-z~^ 
1 I 



V 1 - -2 

(z-1] 



— — 6(1 - z) + 



(5(1 — z) \n.{z„ 



/ln(l - 



w 



(A27) 



11/ -|\ ~\~ '"/ \'~ II LU,dL / I ' 1 -1 

— W\{Z—\)j^ J \ \ — w 

A similar result is obtained for 1/(523^3) and for 1/(523513) except that, in the latter case, 
the singularities occur at w ^ 1 and z ^ Zi. 

Finally, since S23 = sf (1 — w), we make the point that in other cases when the propagator 
l/s23 occurs in the denominator of the matrix elements, it must be combined with the phase 
space factor (1 — and expanded via 



(1-^) 



-l-e 



1 



5{l-w) + 



1 



+ e 



/ln(l - 



w 



0(e) 



(A28) 



The results of this expansion are combined with the phase space factor in Eq. ( |A7| ). To 
ensure in every case that we retain all finite terms and obtain the correct result in the limit 
when /c2 and become coUinear, we always make the full replacement 
sin-2^ 92d92 2^+2^71 tanh^v"' r(l - 2e) 



Jo 



S23 



r2fi 



:i-z)[--jii-w) (--jil-z) 



+ 9{l-z) 

1 

+ 



(l-z) 

( 1 



^ln(^) ^(\n{l-z) 



l-z 



-5 1 -z) + \+6{l-z)[ -\ ^ 



(A29) 
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Most of these terms will vanish since, for example, {1 — z) will usually multiply S{1 — z) 
and (1 — w) will multiply 5(1 — w). In a few special cases, as when there is a Pgg splitting 
function in the coUinear limit, as is the case when the final-state gluon becomes parallel to 
the c-quark in the eg initiated process, the full expansion is needed in order to expose the 
singularity. This singularity may then be factored into the c-quark fragmentation function. 



30 



REFERENCES 



[1] CDF Collaboration, R. Blair et al, FERMILAB-Conf-95/245-E, to be published in 
the Proceedings of the 10th Topical Workshop on Proton- Antiproton CoUider Physics, 
FNAL, May, 1995. 

[2] P. Nason, S. Dawson, and R.K. Ellis, Nucl. Phys. B303, 607 (988); B327, 49 (1989); 
Erratum: B335, 260 (1990); W. Beenakker, H. Kuijf, W.L. van Neerven, and J. Smith, 
Phys. Rev. D40, 54 (1989); W. Beenakker, W.L. van Neerven, R. Meng, G.A. Schuler, 
and J. Smith, Nucl. Phys. B351, 507 (1991); E. L. Berger and R. Meng, Phys. Rev. 
D49 3248, (1994), and references therein. 

[3] P. Aurenche et al, Phys Lett. 140B, 87 (1984); P. Aurenche et al, Nucl. Phys. B297, 
661 (1988); A. P. Contogouris, S. Papadopoulos, and D. Atwood, Theor. Math. Phys. 
87, 374 (1991); H. Baer, J. Ohnemus, and J. F. Owens, Phys. Rev. D42, 61 (1990). 

[4] L. E. Gordon and W. Vogelsang, Phys. Rev. D48, 3136 (1993) and D50, 1901 (1994); 
M. Gliick, L. E. Gordon, E. Reya, and W. Vogelsang, Phys. Rev. Lett. 73, 388 (1994). 

[5] E. L. Berger and J.-W. Qiu, Phys. Lett. B248, 371 (1990); and Phys. Rev. D44, 2002 

(1991). 

[6] P. Aurenche, A. Douiri, R. Baier, M. Fontannaz, and D. Schiff, Z. Phys. C29, 459 

(1985). 

[7] B. Bailey, J.Ohnemus, and J.F.Owens, Phys. Rev. D46, 2018 (1992). 

[8] P. Aurenche, A. Douiri, R. Baier, M. Fontannaz, and D. Schiff, Z. Phys. C24, 309 

(1984). 

[9] E. L. Berger, Phys. Rev. D37, 1810 (1988); M. Mangano, P. Nason, and G. Ridolfi 
Nucl. Phys. B373, 295 (1992). 

[10] S. D. Ellis, Z. Kunszt, and D. E. Soper, Phys. Rev. Lett. 69, 1496 (1992); W. T. Giele, 



31 



E. W. N. Glover, and D. A. Kosower, Phys. Rev. Lett. 73, 2019 (1994); S. D. Ellis and 
D. E. Soper, Phys. Rev. Lett. 74, 5182 (1995). 

[11] E. L. Berger et al, Phys. Rev. D40, 83 (1989); U. Baur et al, Phys. Lett. B318, 544 
(1993). 

[12] E. L. Berger, E. Braaten, and R. D. Field, Nucl. Phys. B239, 52 (1984). 

[13] D.W. Duke and J. F. Owens, Phys. Rev. D26, 1600 (1982); J. F. Owens, Rev. Mod. 
Phys. 59, 465 (1987). 

[14] P. Aurenche et al, Nucl. Phys. B399, 34 (1993). 

[15] M. Gliick, E. Reya, and A. Vogt, Phys. Rev. D48, 116 (1993). 

[16] B. Bailey, E. L. Berger, and L. E. Gordon, Argonne report ANL-HEP-PR-95-87, to be 
published. 

[17] M. Stratmann and W. Vogelsang, Phys. Rev. D52, 1535 (1995). 

[18] R. K. Ellis, M. A. Furman, H. E. Haber, and 1. Hinchliffe, Nucl. Phys. B173, 397 (1980). 

[19] G. AltareUi and G. Parisi, Nucl. Phys. B126, 298 (1977). 

[20] M. Gliick, E. Reya, and A. Vogt, Phys. Rev. D45, 3986 (1992). 

[21] H. L. Lai et al, CTEQ Collaboration, Phys. Rev. D51, 4763 (1995). 

[22] S. J. Brodsky, P. Hoyer, C. Peterson, and N. Sakai, Phys. Lett. B93, 451 (1980); 
S. J. Brodsky C. Peterson, and N. Sakai, Phys. Rev. D23, 2745 (1981). 



32 



Figure Captions 

[1] (a) Lowest order Feynman diagrams for 7 plus c quark production; ki and k2 
are tlie four-vector momenta of the photon and charm quark, (b) Examples of 
virtual corrections to the lowest order diagrams, (c) Examples of next-to-leading 
order three-body final-state diagrams for the gc initial state. 

[2] Charm quark density c{x, Q) as a function of Bjorken x at Q = 10 GeV. The 
solid line shows the expectation of the GRV parton densities [|^, and the dotted 
hue that of the the CTEQ3M densities [0]. 

[3] Cross section da / dpj^dy^ dz as a function of z for p + p— >7 + c + X at ^/s = 1.8 
TeV. We set y"^ = 0. Results are presented in the form of a histogram in bins of 
width A2; = 0.2. In (a), for pj. = 15 GeV, we show the net contribution from the 
lowest order direct process gc 7c and from all the leading order fragmentation 
processes piP2 Psc followed by the collinear fragmentation ^3 7X. In (b) 
and (c), for pj, =15 and 45 GeV, respectively, we display the full cross section 
through next-to-leading order (solid line) and contributions from three important 
0(ag) subprocesses. 

[4] The transverse momentum dependence of da/dpj^dy'^dz, for z integrated over the 
interval 0.2 < z < 2.0. The upper solid line shows the sum of all subprocesses 
through next-to-leading order. The dashed line shows the sum of the 0{as) and 
0(ag) contributions from the eg initial state. The 0(0;^) contributions from the 
gg and cq initial states are shown as dash-dot and dotted curves. The lower solid 
line shows the contribution from the qq (and cc and cc) initial state. 

[5] The renormalization/factorization scale n dependence. For the sum of all con- 
tributing subprocesses, da/dpj'dy'^dz, for y'^ = and z integrated over the inter- 
val 0.2 < z < 2.0, is shown as a function of pj^ for three values of fi/px'- 0.5, 1.0, 
and 2. 
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[6] The renormalization/factorization scale /i dependence of da/dpj^dy^dz. Results 
are shown as a function of 2; at pj. = 20 GeV for three values of n/pr- 0.5, 1.0, 
and 2. 

[7] The K factor defined in the text is shown as a function of for inclusive (i.e., 
non- isolated) photons (solid line) and isolated photons (dashed line); y'^ — and 
0.2<z< 2.0. 
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